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ABSTRACT: The metabolomic profile of the anaerobic fungus Piromyces sp. Fl, isolated from the rumen of goats, and how this is 
affected by the presence of naturally associated methanogens, was analyzed by nuclear magnetic resonance spectroscopy. The major 
metabolites in the fungal monoculture were formate, lactate, ethanol, acetate, succinate, sugars/amino acids and a-ketoglutarate, 
whereas the co-cultures of anaerobic fungi and associated methanogens produced citrate. This is the first report of citrate as a major 
metabolite of anaerobic fungi. Univariate analysis showed that the mean values of formate, lactate, ethanol, citrate, succinate and acetate 
in co-cultures were significantly higher than those in the fungal monoculture, while the mean values of glucose and a-ketoglutarate were 
significantly reduced in co-cultures. Unsupervised principal components analysis revealed separation of metabolite profiles of the fungal 
mono-culture and co-cultures. In conclusion, the novel finding of citrate as one of the major metabolites of anaerobic fungi associated 
with methanogens may suggest a new yet to be identified pathway exists in co-culture. Anaerobic fungal metabolism was shifted by 
associated methanogens, indicating that anaerobic fungi are important providers of substrates for methanogens in the rumen and thus 
play a key role in ruminal methanogenesis. (Key Words: Piromyces sp. Fl, Methanobrevibacter thaueri, Nuclear Magnetic Resonance, 
Citrate) 



INTRODUCTION 

Obligate anaerobic fungi have only been reported in 
gastrointestinal and fecal communities since their first 
discovery in 1975 (Orpin, 1975), and have been found to 
exist exclusively in herbivores, including ruminants and 
hindgut-fermenting animals (Theodorou et al., 1995). These 
fungi play an important role in the primary digestion of 
plant fiber via a broad range of highly active fiber- 
degrading enzymes and their extensive rhizoid system 
(Theodorou et al., 1995). The metabolites of anaerobic 
fungi have so far been reported as formate, acetate, lactate, 
ethanol, succinate, C0 2 and H 2 (Bauchop and Mountfort, 
1981), but no comprehensive metabolomic profiles of 
anaerobic fungi have been reported to date. 

Methanogenic archaea comprise one of the important 
functional groups in the rumen and can utilize hydrogen, 
formate and acetate, metabolites of anaerobic fungi, to 
produce methane (Hungate, 1982). Previous studies showed 
that co-cultures of different anaerobic fungi and 
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methanogens result in increased acetate formation and 
decreased ethanol, lactate, and formate production 
(Bauchop and Mountfort, 1981; Marvin-Sikkema et al., 
1990; Nakashimada et al., 2000), but these co-cultures were 
simple mixtures of axenic pure cultures of an anaerobic 
fungus and a methanogen. In order to provide a more 
natural model system, more closely mimicking the 
interactions of anaerobic fungi and methanogens in the 
rumen ecosystem, Cheng et al. (2006; 2009) obtained a 
mixed natural co -culture of anaerobic fungi and 
methanogens from the goat rumen and found that the 
formate and lactate were reduced in the mixed co -culture 
and acetate was increased. Subsequently, Jin et al. (2011) 
obtained a natural simple co-culture of anaerobic fungus 
and its indigenously associated methanogen from the 
previous mixed co-culture in order to investigate the 
metabolic interaction of anaerobic fungi and methanogens 
in the rumen and their abilities to degrade lignocellulosic 
materials and to produce methane. 

Nuclear magnetic resonance (NMR) spectroscopy is a 
powerful analytical technique for metabolomic studies, 
providing unique and rich information content with high 
sensitivity and reproducibility (Beckonert et al., 2007). It 
has been widely used as a fingerprinting tool for metabolic 
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studies of microorganisms including bacteria and yeast 
(Grivet and Delort, 2009). Metabolites generally associated 
with rumen anaerobic fungi include formate, acetate, lactate, 
succinate, ethanol, C0 2 and H 2 (Theodorou et al., 1996). 
Little information is available on other metabolites of 
rumen anaerobic fungi. In this study, we hypothesized the 
existence of other metabolites produced by rumen anaerobic 
fungi and that the presence of co-cultured methanogens 
could influence the metabolic profiles of rumen anaerobic 
fungi. To test this hypothesis, NMR was used to compare 
the metabolic profiles of three cultures, an anaerobic fungal 
monoculture Piromyces sp. Fl, a simple naturally occurring 
co-culture of Piromyces sp. Fl and its indigenously 
associated methanogen Methanobrevibacter thaueri, and a 
mixed co-culture of anaerobic fungi and methanogens 
derived from rumen digesta. 

MATERIALS AND METHODS 

Microbial cultures 

Three cultures were used in this study: i) a mixed co- 
culture comprising a diverse range of anaerobic fungi and 
methanogens derived from a goat rumen (Cheng et al., 
2006); ii) a natural simple co-culture of an anaerobic fungus, 
Piromyces sp. Fl, and an associated methanogen, 
Methanobrevibacter thaueri, obtained from this mixed co- 
culture using the Hungate roll-tube method (Jin et al., 
2011); iii) a pure culture of Piromyces sp. Fl derived from 
the simple co-culture by adding chloramphenicol to inhibit 
growth of the methanogen. 

Samples collection 

The cultures were grown anaerobically at 39°C in 100 
mL defined medium M2 (Teunissen et al., 1991) 
supplemented with 1% cellobiose. Experimental cultures 
were inoculated by transferring 10 mL of a 3 d-old culture 
to 90 mL medium M2. After incubation at 39°C for 3 d, the 
liquid fractions of the cultures were collected, centrifuged at 
15,500 g and the supernatant were stored at -20°C for NMR 
analysis. 

Chemicals and sample preparation for NMR analysis 

D 2 0 (99.9%) and TSP (3-(trimethylsilyl) propionic- 
2,2,3, 3-c/ 4 -acid, sodium salt) were purchased from 
Cambridge Isotope Laboratories (Miami, FL, USA) and 
Merck (Darmstadt, Germany), respectively. Culture 
supernatant (400 |iL) was mixed with 100 |iL TSP/D 2 0 (1 
mg/mL), centrifuged at 15,500 g for 5 min and the 
supernatant was transferred into NMR tubes. 

NMR analysis 

NMR spectra were obtained on an INOVA 600 MHz 



spectrometer (Varian, USA). The residual water signal was 
suppressed by the pre-saturation method (Lu et al., 2010). 
ID ! H NMR spectra were acquired with the following 
parameters: spectral width 8000 Hz, mix time 0.1 s, 
relaxation delay 2 s, 32K complex data points, 64 transients. 
The data were zero-filled by a factor of 4 and the free 
induction decays (FIDs) were multiplied by an exponential 
weighting function equivalent to a line boarding of 0.5 Hz 
prior to Fourier transformation (FT). The peak of external 
TSP was designated as 0.00 ppm as baseline, and the 
chemical shifts of other chemicals were referenced to the 
peak of external TSP. 

NMR data processing 

The 'H NMR spectra were manually phased and the 
baseline was corrected and referenced to TSP (5 = 0). The 
spectral region (5 = 0.80-8.80) was segmented into 200 
regions at 0.04 ppm per interval. The region from 4.60 to 
5.00 was excluded from the analysis because of the residual 
signal of water in aqueous extracts. The integrated data 
were normalized to the total integrals of each spectrum for 
pattern recognition analysis (Wang et al., 2004) and the data 
were transferred into an excel file. 

Identification of metabolites 

Identification of metabolites was based on their 
chemical shifts and multiplicities with the assistance of 
previous publications and publicly available metabolite 
libraries from the NMR metabolomics database of 
Linkoping, Sweden (MDL) (http://www.liu.se/hu/mdl/ 
main/). 

Data analysis 

*H NMR spectra data were analyzed using the SIMCA- 
P software package (Umetris AB, Umea, Sweden) for 
multivariate statistical analysis. All the variables were 
mean-centered and Pareto-scaled prior to principal 
components analysis. The statistical significance of 
differences between mean values of metabolites was 
analyzed by SPSS software (SPSS Inc., Chicago, USA). 

RESULTS AND DISCUSSION 

NMR spectra and assignment of compounds 

After incubation at 39°C for 3 d, the three cultures grew 
well. Methane was detected in the gas phase of the two co- 
cultures, while no methane was detected from the 
monoculture. 

The NMR spectra of water-soluble metabolites from 
suspensions of three cultures are shown in Figure 1. The 
major metabolites in the anaerobic fungal monoculture were 
formate (5 = 8.46), acetate (5 = 1.94), lactate (5 = 1.34), 
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Figure 1. Nuclear magnetic resonance (NMR) spectra of the water-soluble metabolites from suspensions of anaerobic fungal mono- 
culture (Piromyces sp. Fl; A), a simple co-culture of anaerobic fungus and associated methanogen co-culture (Piromyces sp. 
Fl+Methanobrevibacter thaueri; B) and a mixed co-culture of anaerobic fungi and methanogens derived from rumen digesta, from 
which the simple co-culture was obtained (C). 'H NMR spectra were obtained on an INOVA 600 MHz spectrometer (Varian, USA). The 
'H NMR spectra were manually phased and the baseline was corrected and referenced to TSP (5 = 0). The spectral region (8 = 0.80-8.80) 
was segmented into 200 regions at 0.04 ppm per intervals. The region from 4.60 to 5.00 was excluded from the analysis because of the 
residual signal of water in aqueous extracts. A = Acetate; C = Citrate; E = Ethanol; F = Formate; G = Glucose; L = Lactate; M = 
Unidentified medium components; S = Succinate; K = a-ketoglutarate; S/A = Sugars/amino acids; TSP = Internal reference, 3- 
(trimefhylsilyl) propionic-2,2,3,3-d4-acid. 



ethanol (5 = 1.18), sugars/amino acids (5 = 3.2-4.0) and 
a-ketoglutarate (5 = 4.54) (Kwon et al., 2009). Succinate (5 
= 2.42), a common metabolite of anaerobic fungi, was 
detected in low amounts. Glucose (5 = 5.22), derived by 
anaerobic fungal exoenzymes from cellobiose, was detected 
in large amounts in the culture. 

In the present study, a-ketoglutarate was detected for 
the first time, supporting the assumed metabolic pathway of 
anaerobic fungi proposed by Kwon et al. (2009). Kwon et al. 
(2009) constructed an EST library of rumen fungus 
Neocallimastix frontalis. According to the detected genes 
for glucose metabolism, they elucidated the fungal 
carbohydrate metabolism pathway and postulated that the 



metabolites for Neocallimastix frontalis were acetate, 
ethanol, formate, lactate, succinate and a-ketoglutarate, 
which was oxidized from oxaloacetate by aconitase and 
isocitrate dehydrogenase. In other previous studies, the 
metabolites of anaerobic fungi were reported as formate, 
acetate, lactate, ethanol and succinate. Boxma et al. (2004) 
calculated the relative fluxes of glucose through known 
metabolic pathways and demonstrated that the major 
fermentation end-products of anaerobic fungi were lactate, 
formate, ethanol, acetate and succinate. In the present study 
lactate, formate, ethanol, acetate were detected as major 
fermentation end-products, while succinate was detected in 
low amounts, accounting for 0.16% of the total metabolites. 
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Glucose was detected in the supernatant of the fungal 
monoculture, suggesting that the substrate was degraded 
incompletely in this culture. Large amounts of sugars/amino 
acids were also detected in the supernatant but were 
difficult to separate and identify as reported by Chikayama 
et al. (2008). Chikayama et al. (2008) compiled a 
heteronuclear single quantum coherence-based metabolite 
chemical shift database that contains only NMR spectra of 
standard compounds measured under standardized 
conditions. This database is a compilation of 1,018 'H and 
l3 C chemical shifts corresponding to 142 metabolites, 
showing that many sugars, nucleic acids and organic acids 
were distributed within a narrow range of chemical shifts 
between 3 and 4 ppm for 'H NMR. 

It is interesting that apart from the major metabolites 
formate, acetate, lactate, ethanol and succinate as reported 
previously, citrate (5 = 2.66) was for the first time 
demonstrated as one of the major metabolites in co-cultures. 
In these co-cultures, methanogens can only use simple 
compounds such as H2/CO2, formate, acetate etc. to produce 
methane, and Methanobrevibacter sp. strains extensively 
use H 2 /C0 2 or/and formate as substrate to produce methane 
(Balch et al., 1979). But the mechanism of citrate 
production by anaerobic fungi is unknown and needs to be 
further investigated. Meanwhile, two peaks at 4.14 ppm and 
3.66 ppm chemical shift, representing lactate and ethanol 
respectively, were observed in the co-cultures. Other peaks 
at 7.38 ppm, 4.42 ppm, 2.78 ppm, 1.58 ppm and 1.26 ppm, 
were not observed in the fungal monoculture, and were 
found, by NMR analysis, to come from the growth medium 
(data not shown). 

The mean values of the metabolites from three co- 
cultures are shown in Table 1; suggesting that the 
concentrations of formate, acetate, lactate, citrate, succinate 
and ethanol were significantly higher in co -cultures than the 



anaerobic monoculture. But, there was no statistical 
difference between the concentrations of these metabolites 
in the two co-cultures, implying similar fungal metabolism 
in the two co-cultures. In agreement with the suggestion by 
Kwon et al. (2009), a-ketoglutarate was one of the major 
metabolites in the fungal mono-culture. However, the 
concentration of a-ketoglutarate was similarly low in the 
two co-cultures, significantly lower than in the anaerobic 
fungal monoculture. Meanwhile, sugars/amino acids within 
3.0 to 4.0 ppm, except ethanol (8 = 3.66 ppm), were found 
at much lower concentrations in the co -cultures than in the 
fungal monoculture. Glucose split by anaerobic fungal 
exoenzymes from cellobiose was also detected at 
significantly lower concentration in co-cultures than the 
fungal monoculture. This may imply that the substrate for 
anaerobic fungi was utilized completely in the co-cultures. 
Consequently, the production of metabolites such as 
formate, acetate, ethanol, lactate and succinate was 
increased in these co-cultures. 

Bauchop and Mountfort (1981) reported that the amount 
of acetate in culture supernatant was much higher when the 
anaerobic fungus Neocallimastix strain PN1 was co- 
cultured with Methanobrevibacter sp. strain RA1 than in the 
fungal monoculture. The authors subsequently reported that 
when Methanosarcina barkeri was cultured with 
Neocallimastix strain PN1 and Methanobrevibacter sp. 
strain RA1 co-culture, acetate could not be detected in the 
tri-culture (Mountfort et al., 1982), suggesting that the 
amount of acetate produced in the co -culture of anaerobic 
fungi and methanogens was dependent on the co-cultured 
methanogen species. However, these authors reported that 
lactate, ethanol and formate were produced at much lower 
levels in the co-culture than in fungal mono-culture 
(Bauchop and Mountfort, 1981; Mountfort et al., 1982). 



Table 1. Mean values of detected metabolites in an anaerobic fungal mono-culture and two fungal/methanogen co-cultures 



Chemical shifts (5) 1 
(ppm) 



Compounds 



Cultures 



B 



1.18 


Ethanol 


291.5±20.3 a>4 


489.2±24.4 b 


429.9±51.6 b 


1.34 


Lactate 


186.2±8.0 a 


1984.9±124.5 b 


1434.9±436.7 b 


1.94 


Acetate 


204.2±10.2 a 


710.4±2.6 b 


653.5±94.8 b 


2.42 


Succinate 


15.8±l.l a 


75.7±0.7 b 


65.7±11.4 b 


2.66 


Citrate 


19.2±1.9 a 


283.0±18.2 b 


248.0±10.1 b 


4.54 


a-Ketoglutarate 


145.4±5.4 a 


27.1±0.3 b 


39.2±10.4 b 


5.22 


Glucose 


112.5±5.4 a 


9.8±0.5 b 


23.3±10.2 b 


8.46 


Formate 


115.6±2.1 a 


429.3±16.1 b 


315.7±62.7 b 



The baseline was corrected and referenced to TSP (3-(trimethylsilyl) propionic-2,2,3,3-rf 4 -acid, sodium salt; 5 = 0). 

2 Identification of metabolites was based on their chemical shifts and multiplicities with the assistance of previous publications and publicly available 
metabolite libraries from the NMR metabolomics database of Linkoping, Sweden (MDL) (http://www.liu.se/hu/mdl/main/). 

3 A = Anaerobic fungal monoculture; B = Co-culture of anaerobic fungus and methanogen; C = Mixed co-culture of anaerobic fungi and methanogens 
derived from the rumen. Mean values are presented as mean+SEM. 

4 The statistical significance of differences between mean values of metabolites was analyzed by SPSS software (SPSS Inc., Chicago, USA). In each row, 
values followed by different superscripts are statistically significant (p<0.05). 
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Figure 2. PLS-DA and O-PLS analysis derived score plots of the first two principal components (PCI and PC2) of suspensions of an 
anaerobic fungal mono-culture (al-3), a simple co-culture of an anaerobic fungus and methanogen (bl-3) and a mixed co-culture of 
anaerobic fungi and methanogens derived from goat rumen digesta (cl-3). The 'H NMR spectra data were analyzed using the SIMCA-P 
software package (Umetris AB, Umea, Sweden). All the variables were mean-centered and Pareto-scaled prior to principal components 
analysis and the ellipse represents the 95% confidence region. 



Similarly, Joblin et al. (1990) reported that levels of lactate 
and formate were much lower in a co-culture of N. frontalis 
and Methanobrevibacter smithii than fungal mono-culture. 
Cheng et al. (2009) also reported a much lower lactate and 
formate concentrations in mixed co -culture of anaerobic 
fungi and methanogens than fungal culture. These reports 
contrast with data reported here, which demonstrated that 
formate, acetate, ethanol and lactate all increased in co- 
cultures. These differences may be due to the different 
substrates used in different studies. Jin et al. (2011) 
investigated the same co-cultures as that used in this study 
and reported that there were no significant differences on 
the levels of lactate and ethanol between co -culture and 
fungal culture, when growing on corncob. Our other study 
on the same co-culture showed that the acetate, ethanol and 
lactate all significantly increased in co-cultures, when 
growing on cassave (Liu, 2009). Thus, this metabolic 
differences on different substrates might be dependent on 
anaerobic fungal and methanogen species. 

It should also be noted that in previous studies, culture 
metabolites were assayed by chemical methods that 



analyzed individual compounds (Mountfort et al., 1982; 
Cheng et al., 2009), while the NMR analysis used here is 
more sensitive and reliable. NMR has been used to 
investigate microbial physiology and metabolism for over 
30 yrs (Grivet, 2001). Carrieri et al. (2009) used this 
approach to demonstrate that the water-soluble metabolites 
of two cyanobacterial species Arthrospora maxima and 
Synechococcus 7002 were formate, lactate, ethanol, acetate 
and succinate. Pham et al. (2005) identified the secondary 
metabolites from Streptomyces violaceoruber and results 
showed that more than 50% of the identified metabolites 
were new compounds. 

Principle components analysis 

The data were analyzed with PLS-DA and O-PLS 
(Figure 2A and 2B), but only O-PLS could discriminate the 
three cultures, and the first two principal components 
cumulatively accounted for 86.3% of the total variation. 
O-PLS was then used to investigate the major compounds 
that contributed to discrimination, and an O-PLS derived 
variable importance plot is shown in Figure 3. The most 




Variable ID (chemical shift) 

Figure 3. O-PLS derived variable importance plot of the suspensions of an anaerobic mono-culture and two fungal/methanogen co- 
cultures. The variable ID corresponds to a chemical shift (ppm). 
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important compound contributing to the discrimination was 
lactate, corresponding to a chemical shift 1.34 ppm, 
followed by acetate, sugars/ AA, ethanol and formate (Table 
2). 

Metabolic pathway of anaerobic fungi with/without 
methanogens 

According to previous reports, lactate, formate, ethanol 
and succinate were formed in the cytosol of anaerobic fungi, 
and acetate, formate and H 2 /C0 2 were formed in the 
hydrogenosome (O'Fallon et al., 1991; Akhmanova et al., 
1999; Boxma et al., 2004). In a recent report by Kwon et al. 
(2009), a-ketoglutarate was proposed to be produced in the 
cytosol. In the present study, citrate was reported as 
metabolite of anaerobic fungi in co -culture with 
methanogens. Based on these reports, a modified metabolic 
pathway of glucose by anaerobic fungal monoculture is 

Glucose 



Table 2. The five major compounds contributing to O-PLS 
clustering of an anaerobic fungal mono-culture and two 
fungal/methanogen co-cultures 



No. 


Chemical shift (5) 
(ppm) 


VIP 1 value 


Corresponding 
metabolites 


1 


1.34 


5.28 


Lactate 


2 


1.94 


3.11 


Acetate 


3 


3.50 


3.09 


Sugars/Amino acids 


4 


3.94 


3.09 


Sugars/Amino acids 


5 


3.42 


2.79 


Sugars/Amino acids 



1 VIP = Variable importance plot. 

shown in Figure 4. 

Anaerobic fungi possess both reductive and oxidative 
branches of the TCA cycle in the cytosol, producing 
succinate and a-ketoglutarate respectively as proposed by 
Kwon et al. (2009). In the present study, succinate was 



Glucose 
Glucose-6-Pi 



|GPI 

Fructose-6-Pi 
^PFK 

Fructose- 1,6-Pi 

Jfba 

Glyceraldehyde-3-Pi 
^GAPDH 
1 ,3-bisphosphoglycerate 
|PGK 
3 -Phosphoglycerate 
^IPGM 

2-Phosphoglycerate 

^Enolase 

Phosphoenol pyruvate 

I PK 
LDH ■f 
Lactate ^ Pyruvate 

Formate ^ 



iPFL 
_ -CoA 



Acetyl-CoA 

Acetalcfehyde 

^ADHE 
Ethanol 



Cytosol 



Acetate 

scsb| 

Acetyl-CoA 



Hydrogenosome 



-► Formate 



Pyruvate 

me| 

Malate 



-+- C0 2 , H 2 



PEPCK 



>. Oxaloacetate MDH » Malate ■ 



_^ Fumarate 



FR 



^. Succinate 



Reductive branch 



Citrate 



AC 



Aconitate 



AC ICD 
^. Isocitrate a-Ketoglutarate 

Oxidative branch 



Figure 4. Metabolic pathway of carbohydrate by anaerobic fungi. Metabolic pathway of carbohydrate by anaerobic fungi based on 
previous reports by O'Fallon et al. (1991), Akhmanova et al. (1999), Boxma et al. (2004) and Kwon et al. (2009). AC = Aconitase; 
ADHE = Aldehyde/alcohol dehydrogenase; F = Fumarase; FBA = Fructose-bisphosphate aldolase; FR = Fumarate reductase; GAPDH = 
Glyceraldehydes 3-phosphate dehydrogenase; GPI = Glucose-6-phosphate isomerase; HK = Hexokinase; ICD = Isocitrate 
dehydrogenase; IPGM = 2,3-bisphosphoglycerate-independent phosphoglycerate mutase; LDH = Lactate dehydrogenase; MDH = Malate 
dehydrogenase; ME = Malic enzyme; PEPCK = Phosphoenolpyruvate carboxykinase; PFK = Phosphofructokinase; PFL = Pyruvate 
formate lyase; PGK = 3-phosphoglycerate kinase; PK = Pyruvate kinase; SCSB = Succinyl-CoA synthetase, beta subunit. 
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detected in lower amounts than a-ketoglutarate in anaerobic 
fungal monoculture, suggesting that the reductive branch 
was less active than the oxidative branch of the TCA cycle 
in the cytosol of fungal monoculture. However, succinate 
production increased and a-ketoglutarate decreased 
dramatically in the co-cultures, implying that the reductive 
branch was stimulated, and the oxidative branch was 
inhibited, resulting in the accumulation of citrate. The 
citrate synthase, converting oxaloacetate to citrate, has not 
been detected until now, even though Kwon et al. (2009) 
postulated that there might be a possible existence of citrate 
synthase. The co-culture system in this study might be a 
useful tool for investigating the citrate formation in 
anaerobic fungi. Meanwhile, the oxidation reactions in the 
cytosol of anaerobic fungi and effects of co-cultured 
methanogens on these reactions need much more 
investigation. 

In the hydrogenosome of anaerobic fungi, the metabolic 
pathway was also stimulated as revealed by increased 
acetate production. Results in this study showed that the 
acetate production in the co-cultures was twice that in the 
monoculture. According to previous reports, the H 2 formed 
in the hydrogenosome was provided for the growth of co- 
cultured methanogens, and cannot be detected in the co- 
cultures (Marvin-Sikkema et al., 1990; Teunissen et al., 
1992). Until now, the knowledge of the metabolic pathway 
in the hydrogenosome of anaerobic fungi was limited. And 
thus it is difficult to analyze the effects of associated 
methanogens on the metabolic pathway in the 
hydrogenosome of the fungi. In the future, the metabolism 
in the hydrogenosome of anaerobic fungi needs to be 
further investigated in mono- and co-cultures. 

This study was the first to investigate the metabolism of 
anaerobic fungi by metabolomics and the novel finding of 
citrate as one of the major metabolites of anaerobic fungi 
co-cultured with methanogens by NMR may suggest a new 
yet unknown pathway in the co-culture. However, the 
metabolic pathway of anaerobic fungi, especially in the 
presence of methanogens, needs to be further investigated 
in the cytosol as well as in the hydrogenosome. The present 
study also showed that the anaerobic fungi were important 
electron contributors for methanogens in the rumen and 
their contribution to ruminal methane emission needs to be 
assessed. 
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